Robust topological insulator surface conduction under strong surface disorder 
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Topological insulators are characterized by specially protected conduction on their outer bound- 
aries, i.e. on their surlaces. We show that the protected surface conduction exhibited by 2-D 
topological insulators is independent of non-magnetic surface disorder. In particular, the surface 
band remains conducting even when surface state inhomogeneities destroy the characteristic linear 
Dirac relation between energy and momentum. The main effects of disorder are to pull the surface 
states into the disordered layers, decrease their Fermi velocity, and increase the density of states. 
These effects are controlled by a resonance between the disorder potential and the bulk bands. The 
resonance's energy is set by the bulk band width; protection of the Dirac cone is controlled by the 
bulk band width, not the bulk band gap. 
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I. INTRODUCTION 

Recently a new kind of material has been predicted 
and measured: topological insulators, which do not per- 
mit current to flow through their interior but do al- 
low metallic conduction along their surfaces .1^^ The con- 
ducting surface states are protected topologically, mean- 
ing that they are safeguarded by the bulk's insulating 
property from any finite-size perturbations. The surface 
states obey a linear Dirac dispersion like that seen in 
graphene, but unlike graphene their spin degeneracy is 
strongly broken, the angles of the spin and the momen- 
tum are locked to each other, and T symmetry prohibits 
any gap in the Dirac dispersion. Topological insulators 
(TIs) have won tremendous theoretical and experimental 
interest because these special properties have potential 
for spintronics and quantum computation. However TIs 
continue to face critical experimental challenges. Chief 
among these are elimination of residual bulk carriers from 
dopants introduced by the TI growth process and by ex- 
posure to air, control of the Fermi level which responds 
strongly to doping, and fabrication of clean and stable 
surfaces. 

This present article focuses on how non-magnetic dis- 
order located only on the TI surface affects topological 
protection. Several experiments have shown that TI sur- 
face states are not destroyed even under exposure to 
aij-IMIl However three recent experiments measured the 
surface band while progressively doping the TI surface 
with impurities .l^^li^ They saw the angularly resolved 
photoemission spectroscopy (ARPES) signal - an image 
of the Dirac cone - become progressively fainter and more 
blurred as impurities were added. Moreover several re- 
cent experiments have modified the TI surface by in tro- 
ducing an Al capping layer, or by gating the TI.mAn- 
other experiment reported that exposure to air causes 
growth of an oxide layer on top of the Tl^^ Therefore a 
better understanding is needed of surface disorder's ef- 
fects on topological protection of TI surface states. In 
particular, we will see that direct measurements of the 
Dirac cone (such as ARPES and Schlubnikov-de Haas 



measurements of 3-D TIs) are quite sensitive to disorder. 
Strong surface disorder can disrupt the Dirac cone and 
extinguish these experimental signals while leaving the 
surface conductance intact. 

The current literature offers two lenses for understand- 
ing topological protection. The first lens is focused on 
bulk properties which protect the surface states. In TIs 
the bulk valence and conduction bands are split both by a 
strong spin-orbit interaction and by a gap whose width is 
controlled by a mass parameter m. If the mass and spin- 
orbit terms have opposite signs, the conflict between the 
two signs entangles the valence band with the conduc- 
tion band, and the states of each band acquire a special 
texture centered on the fc = F point.'^^ This situation, 
called band inversion, strongly affects the TI's electronic 
structure both at its surfaces and also at its boundaries 
with any material having a positive (non-inverted) mass. 
On these surfaces the band inversion creates a band of 
surface states which connects the valence and conduction 
bands with each other and therefore bridges the bulk gap. 
If the Fermi level lies in the bulk gap, then any states in 
the bulk decay exponentially and current flows only in 
the surface states caused by the band inversion. 

The surface states are safeguarded - topologically pro- 
tected - by the bulk band gap. Destruction of the sur- 
face band can occur only when surface states tunnel 
from one side of the sample to the other, i.e. when 
the bulk states at the Fermi level become delocalized. 
Such del ocaliz ation could be caused by the sample's 
geometryj ^*' ! ^^ ! by interaction s,!^ o r by disorder-assisted 
percolation through the bulk.l^fllSil Whatever its source, 
bulk delocalization is not possible unless the bulk gap 
is destroyed throughout the sample, which is not pos- 
sible with surface disorder. In particular, surface dis- 
order leaves the bulk band gap unchanged. Therefore 
we should expect TI surface states to be robust against 
surface disorder. If the surface disorder is strong enough 
to locally change the band structure in the disordered re- 
gion, then the TI surface states will simply detour around 
the disorder and ride the boundar y betw een the disor- 
dered surface and the clean TI buIk.l^SlSS] 
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The second lens for understanding topological protec- 
tion focuses on the surface states and their spin texture. 
T symmetry requires that near the T point the surface 
band of a topological insulator must be a linear Dirac 
cone, and that the angles of the spin and the momentum 
must be locked together. This spin texture has far reach- 
ing consequences for surface conduction both in two and 
in three dimensions. In a 2-D TI the Dirac cone is com- 
posed of only two edge states which have opposite spin 
quantum number and opposite momenta |p| = zLEp/vp, 
where Ep and vp are respectively the Fermi energy and 
Fermi velocity. Non-magnetic disorder is unable to flip 
the spin and therefore leaves the edge conductance always 
quantized at 2Go, where Go is the conductance quantum. 
In a 3-D TI the spin's effect is bit more subtle, since it 
permits scattering at any angle except 180 degrees. In 
this case spin-momentum locking implies that electrons 
pick up a TT Berry phase every time that they make a 
closed loop. This Berry phase reverses the sign of dis- 
order's effects on the conductance: weak disorder causes 
an increase in the conductance, as is typical in systems 
with strong spin-orbit coupling.^'* Non-perturbative nu- 
merical studies of purely surface conduction with a single 
Dirac cone surface band indicate that the conductance in- 
creases wit hout l imit, proportional to InL where L is the 
wire length.l^SISSl This is much different from long non- 
topological wires, where the conductance invariably de- 
creases to zero. These arguments based only on the Dirac 
cone and on spin-momentum locking give a detailed pic- 
ture of how topological protection affects scattering pro- 
cesses. 

How does topological protection work when surface 
disorder becomes large? Recently Schubert et aPSl per- 
formed a numerical simulation of 3-D TIs with surface 
disorder and found that the surface band evolves through 
three regimes when the surface disorder is increased. At 
small disorder the surface states are close to plane waves 
and follow a linear Dirac dispersion. However when the 
disorder grows larger the in-gap states are no longer ho- 
mogenous plane waves and therefore the linear Dirac re- 
lation between momentum and energy disappears. This 
is consistent with several experimental studies which saw 
the Dirac cone weaken and blur as surface disorder was 
increased. E^^^ Lastly at very large surface disorder the 
in-gap states move just inside the disordered region and 
return to the Dirac dispersion found at weak disorder. 
Based on these numerical results, Schubert et al argued 
that when the disorder is larger than the gap the bulk and 
surface states mix, the Dirac cone is destroyed, and the 
TI surface states are no longer topologically protected. 

The present paper studies the same physics, but in 2- 
D TIs, where the quantized 2Gq conductance simplifies 
analysis, much larger sample sizes are accessible, and sta- 
tistical errors and finite size effects are more easily con- 
trolled. We confirm that there are three disorder regimes. 
Like Schubert et al, we show that in the intermediate dis- 
order regime the surface states are inhomogenous (not 
plane waves) and therefore do not follow any dispersion 



relation between energy and momentum. However we 
also show that surface conduction is topologically pro- 
tected at any disorder strength - the 2Gq quantization 
is never destroyed by surface disorder. In this interme- 
diate disorder regime direct measurements of the energy 
dispersion will find no signal of the Dirac cone, while 
transport experiments will show robust surface conduc- 
tion. 

Going further, we find that the inhomogeneous states 
at intermediate disorder are controlled by a resonance 
between the disorder potential and the bulk bands. Both 
the resonance and the accompanying disruption of the 
Dirac cone occur at the characteristic energy scale of the 
bulk bands; the Dirac cone is protected by the bulk band 
width, not by the much smaller bulk band gap. The main 
effects of disorder in the resonant region are to pull the 
surface states into the disordered outer layers, increase 
the density of states, and decrease the Fermi velocity. 

The structure of this paper is as follows. In section 
O we describe our model of a 2-D TI and the numeri- 
cal methods which we use to study it, in section [HT] we 
present and discuss our numerical results, and in section 
IV we discuss the interpretation and physics of these re- 
sults. 



II. METHODS 

We study the Bernevig-Hughes-Zhang (BHZ) tight 
binding model^ of 2-D topological insulators. This 
model was obtained by starting with the six-band Kane 
model of the Tq and Fs bands and then removing the LI 
sub-band which splits off. The remaining basis has four 
orbitals: the first two are \s,Sz = 1/2), \p,Sz — 1/2), 
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and the last two have Sz reversed 
momentum representation is: 



The BHZ model's 
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512meVnm^, B — 686meVnm^, 
364:.5meVnm, a — bnm, M = —IQmeV (1) 



^x,y,z are the Pauli matrices, and the sines and cosines 
in the Hamiltonian ensure the correct topological band 
structure, a = 5 nm is the lattice spacing. The penetra- 
tion depth of the in-gap surface states is always less than 
5a 25 nm, which gives us good control of finite size ef- 
fects when calculating the conductance in strips of width 
3a to 28a. Finite size effects are practically absent in our 
other calculations, where we have checked that changing 
the system size from 200a x 200a to 400a x 400a makes 
very little difference in our results. The mass M can be 
tuned by changing the sample thickness, and we choose 
a negative mass M = — 10 meV so that there is a band 
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FIG. 1: (Color online.) Effects of disorder on the spectral density A{k,E) and typical surface states. Panes a, b, and c show 
the spectral density as a function of both energy E and wave- number k at small disorder (W = 20 meV), intermediate disorder 
{W = 100 meV), and large disorder {W — 2000 meV). The Dirac cone's two lines are clearly visible at small disorder and large 
disorder, but disappear in a range starting at Vl^ = 80 meV and extending to near W = 2000 meV. Pane d shows the averaged 
spectral density on a cross section E = —2.5 meV; the Dirac cone is signaled by two symmetric peaks and is absent in a range 
starting at = 80 meV. Panes e-h show typical surface states at four representative disorder strengths. In all cases they are 
extended and conducting, but a.t W — 100, 1000 meV they lose their plane wave character and therefore are not visible in the 
spectral density. 



inversion. This mass can produced by fabricating a quan- 
tum well with a thickness of ^ 70 A.l^ With this mass the 
bulk exhibits a band gap in the interval [—10, 10] meV, 
and the Dirac point lies near 8 meV. The valence band 
extends over [-62,-10] meV, and the conduction band 
extends over [10, 373] meV. The bulk bandwidths of 52 
meV, 363 meV are important - we will show that they set 
the scale of the disorder's effects on the Dirac dispersion. 

The BHZ model conserves the Sz component of the 
spin, which is protected by bulk inversion symmetry. 
This symmetry can be broken by a quantum well or gate 
electrod^^, and also can be broken by disorder. We in- 
troduce disorder on the outer two layers (2a = lOnm) 
of the TI, and this disorder is white noise without a ny 
correlation between sites. We follow several papergSSEil 
which preserved the Sz symmetry and used nonmag- 
netic on-site disorder randomly distributed in the interval 
[—W/2, W/2], where W is the disorder strength. With 
this disorder the BHZ model factorizes into two indepen- 



dent 2x2 Hamiltonians which are time-reversal coun- 
terparts of each other, and its universality class changes 
from symplectic to unitary. Althou gh t his change is per- 
tinent to disorder effects in the bullpi"2S!, it does not af- 
fect the present paper's results on surface disorder with a 
clean bulk. Here the conducting states in the gap are edge 
states controlled a 2 x 2 effective Hamiltonian whose sym- 
metry class is independent of the BHZ model's 5*2 sym- 
metry. All of our calculations of the conductance, eigen- 
vectors, and eigenvalues concern themselves only with 
the Sz — 1/2 sector, which is a time-reversed copy of the 
5^ = — 1/2 sector. 

With the exception of our conductance data, all of 
our numerical results come from a leads-free geometry. 
However for the conductance we use a strip with two at- 
tached clean semi-infinite leads that have widths equal 
to the width of the strip itself. We evaluate the conduc- 
tance using the Caroli formuliPSl^Zl g = -'^Tr{{T,l - 
I]i)G£^(E^ - E^JG^x), where G", C are the advanced 
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and retarded Green's functions connecting the left and 
right leads and Si,_R are the self-energies of the leads. 
We evaluate the lead self-energies using the well-known 
iterative technique developed by Lopez Sancho et alP^ 

III. RESULTS 

Panes a-d of Figure [T] summarize the spectral den- 
sity - the square A{k,E) = - of 
the eigenstates {ipn), resolved in momentum space. This 
quantity is very attractive because in pure samples it 
shows a clear X-shaped signal marking the surface state 
Dirac cone, and because in 3-D TIs it may be measured 
with ARPES. Such measurements have been very popu- 
lar both for proving that particular materials are topo- 
logical insulators and for diagnostic evaluations of in- 
dividual experimental samples. In panes a-c of Fig. [l] 
we report the spectral density in the band gap of a sin- 
gle 2-D 400a X 400a sample, using a logarithmic color 
scale. Pane a shows the spectral density of a weakly dis- 
ordered sample {W = 20 meV), which is substantially 
the same as the density of a disorder-free sample. The 
two straight lines of the Dirac cone are very clear, and 
the Dirac point lies near 8 meV. Above 10 meV and be- 
low — 10 meV the bulk valence and conduction bands 
are visible as very bright curves. In pane b we increase 
the disorder strength to W = 100 meV and observe that 
the Dirac dispersion is no longer visible while the bulk 
bands remain unchanged. Lastly pane c shows that when 
the disorder strength is far larger than the band width 
{W = 2000 meV vs. 363 meV) the surface states again 
manifest themselves in a clearly visible linear Dirac dis- 
persion. Pane d tells the same story with more precision 
by focusing on a cross-section of A{k, E) at E = —2.5 
meV and averaging over 1000 200a x 200a samples. The 
data for small disorder {W = 20, 70 meV) show two peaks 
at |fc| = \E — Eol/hvp corresponding to the two branches 
of the Dirac cone. The peaks disappear at = 80 meV 
and reappear near W = 2000 meV. These results agree 
completely with Schubert et al's study of the spectral 
density in 3-D TIs^'^ which found that the Dirac cone is 
absent from the spectral density at intermediate disor- 
ders but reappears at large disorders. Our results also 
are consistent with several experimental ARPES obser- 
vations of a weakened Dirac cone in 3-D TIs with surface 
disorder .ESHU 

To understand the fate of the surface states we calcu- 
late individual surface states at four representative dis- 
order strengths and show their probability distributions 
[■(/'(a?)!^ in panes e — /i of Figure [l] We plot the lower 
5 layers of 200a x 200a samples, choose E = 7.5 meV, 
and employ a linear color scale. In all cases the great 
majority of surface states remain extended along the en- 
tire perimeter of the TI sample and strongly avoid the 
sample's interior. By collecting large numbers of surface 
states we have found a small sub-population of very lo- 
calized states in the gap, but these states are a small 
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FIG. 2: (Color online.) Robust topological protection of sur- 
face conduction, limited only by tunneling through the bulk. 
Pane a shows how we determine the tunneling length Lt- The 
average conductance follows a universal curve, which is shown 
in the left-most curve. The thin black line shows a best fit 
with the hyperbolic tangent function. We determine the tun- 
neling length Lt by finding the point where the average con- 
ductance is equal to 1. Pane b shows Lt as a function of strip 
width w for several Fermi energies and disorder strengths, 
including W = 10* meV which far exceeds the band width. 
Each curve converges rapidly to a straight line, proving that 
the tunneling length becomes exponentially large when the 
strip is widened, and that surface conduction is topologically 
protected in wide strips. 

fraction of the total."^ Pane e shows a surface state at 
zero disorder. It is a plane wave; its density is constant 
along the sample edge. Likewise pane h shows that at 
large disorder {W — 10^ meV) the surface states shift in- 
ward into the sample, avoiding the outer two disordered 
layers, and therefore have a smooth plane wave charac- 
ter. However at intermediate disorders W = 100, 1000 
meV the surface states are strongly affected by the disor- 
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FIG. 3: (Color online.) The resonance between disorder and the bulk band which is responsible for disrupting the Dirac cone. 
Pane a shows the density of states at seven Fermi levels inside the band gap. A clear resonance is visible at the characteristic 
bulk band widths [52, 363] meV, which are marked with vertical dashed lines. The green dashed line plots the profile of the 
Fermi velocity at i? = meV, which is inversely proportional to the DOS. The units are arbitrary. Pane b shows the average 
distance from the edge (d), while pane c plots the probability as a function of distance from the edge. Both panes show 

that in the resonant region the surface states are pulled into the disordered layers. The inset in pane c shows that deep inside 
the bulk the surface states are controlled by the penetration length plus a shift caused by the resonance. Pane d shows the 
density of states when the band gap lies in the interval [—10, 10] meV, or instead in the interval [—5, 5] meV. The curve far to 
the left shows where the [—5, 5] meV DOS curve would be if the resonance were controlled by the band gap instead of the band 
width. 



der and become very inhomogeneous. This explains the 
absence of a Dirac cone structure in the spectral den- 
sity. Extended surface states are present regardless of 
the disorder strength, but at = 100, 1000 meV they 
are not plane waves and therefore have no structure in 
momentum space. 

Although the surface states are no longer visible in 
the spectral density, they remain topologically protected 
and in wide strips the conductance remains forever quan- 
tized. In Figure [2] we calculate the conductance in long 
TI strips and show that the only mechanism which de- 
stroys conduction is tunneling across the strips. In these 
conductance calculations we introduce disorder only in 
the outermost layer of the TI. The tunneling occurs at 
a characteristic strip length Lt which can be measured 
by calculating the average conductance as a function of 
strip length L, as seen in Fig. [2^. We calculated this 
function at a variety of disorder strengths and Fermi 
energies, averaging over 3000 disorder realizations and 
keeping the strip width fixed at = 20a. In each case 
(except for small strip widths) the average conductance 
follows an identical curve, as shown in the left-most curve 
where we have shifted all the data to coincide with each 
other. The universal curve is nearly linear in the interval 
G = [0.4, 1.4] and is fairly similar to the hyperbolic tan- 
gent, which is marked by a thin black line. This makes 
us confident that we can reliably calculate the tunnel- 
ing length scale Lt by determining the value of L which 
satisfies G{L) = 1. In pane b we plot Lf as a function 
of strip width w. Each curve is determined by varying 
the strip width while keeping the disorder strength and 
Fermi level constant. We explore the bulk band gap be- 
tween [—10, 10] meV by calculating four Fermi energies 
E = -7.5,-5,0,5 meV. We also explore three disorder 



strengths W — 10^, 10'^, 10* which all far exceed the bulk 
band gap, and the last of these far exceeds the band 
width. At small strip widths the curves show some bend- 
ing caused by finite size effects. Nonetheless each curve 
quickly converges to a straight line; the tunneling length 
Lt is proportional to Lt oc e^' and is exponentially large 
in a wide strip. This proves that the conductance remains 
forever quantized in an infinitely wide strip; surface con- 
duction is robustly protected regardless of the disorder 
strength. The surface states inside the bulk band gap 
remain always extended and conducting. 

Lastly we study the physics which disrupts the Dirac 
cone signal in the spectral density, which is a resonance 
between the disorder potential and the bulk bands at 
the characteristic bulk band widths [52, 363] meV. Pane 
a of Figure [3] reports the density of states (DOS)'s de- 
pendence on disorder, which we obtained by collecting 
eigenvalue-eigenstate pairs near specified Fermi energies 
E — —7.5, ...,7.5 meV and averaging over 100 disorder 
realizations. First we note that the DOS (and all other 
observables) evolves continuously; the absent Dirac cone 
can not be explained by an absence of surface states. In- 
side the interval E — [52, 363] meV the DOS first grows 
sharply and then begins to decay more slowly. The peak 
over base ratio of this resonance is about 5.8. On the left 
shoulder of the resonance we see a spread between the 
DOS curves measured at different Fermi energies, indi- 
cating that the density of states increases first near the 
valence band and soon after at higher Fermi energies. 
The valence band has a much smaller bandwidth, so it 
is natural that the left shoulder of the resonance has a 
predominantly valence band character. Outside of the 
resonance the DOS tends toward its disorder-free value. 
Except for in the resonance's left shoulder which is rela- 



6 



tively small, the DOS is independent of the Fermi energy, 
as expected of Dirac states on the edge of a 2-D TI. 

Panes b and c of Figure [3] show that the resonance pulls 
the edge states into the disordered sites on the bound- 
ary of the TI. Two authors have a lready briefly discussed 
this behavior at a qualitative leveP^'^ here we measure 
it quantitatively and establish its dependence on the dis- 
order strength and the Fermi energy. Pane b plots the 
average distance from the edge - the expectation value 
(d) = / (Px\'ijj{x)\'^d{x) of the distance d{x) to the near- 
est edge. Again we have averaged over 100 disorder real- 
izations. Outside of the resonant interval [52, 363] meV 
(d) depends strongly on the Fermi level and is maximized 
near the Dirac point E = 8 mcV. This behavior is linked 
to the penetration depth which is also energy dependent, 
and contrasts strongly with the resonant interval where 
(d) is practically independent of energy. Inside the reso- 
nance the state's position is determined by disorder and 
not by the penetration depth, (d) is strongly reduced to 
about one lattice unit; the surface states reside almost en- 
tirely in the outer two layers where we have introduced 
disorder. We see this more precisely in pane c, which 
plots the probability P{d) = J d'^x\^{x)\'^ S{d—d{x)) that 
the state is located on layer d. We examine four disor- 
der strengths, average over 10^ disorder realizations, and 
keep the Fermi energy fixed at 5 meV. Inside the reso- 
nant interval P{d) becomes strongly concentrated on the 
outer two layers; the probability on the outer layer P(0) 
almost doubles from the pure case to = 100 meV. As 
the disorder is increased beyond the band width the sur- 
face states are expelled from the disordered layers which 
are no longer part of the TI because their band structure 
is overwhelmed by disorder. The surface states concen- 
trate at the boundary between the disordered layers and 
the clean layers which retain their TI band structure. 
The inset shows that deep inside the TI P{d) is always 
controlled by the penetration depth and decays exponen- 
tially. Here the only effect of disorder is to shift P{d) first 
inward at intermediate disorder strengths and then out- 
ward at large disorder. This effect is similar to the phase 
shifts found in the theory of scattering events. 

The above discussion suggests that the bulk band 
width not the bulk band gap controls the energy scales 
at which Dirac cone is protected or destroyed. This re- 
sult contrasts with at least one previous papeP^ which 
expected that the Dirac cone is protected by the bulk 
band gap. Therefore we checked this question directly by 
changing the gap to [—5, 5] meV, in contrast to the rest of 
our calculations where we set the gap to [—10, 10] meV. 
We recomputed the DOS, the average distance from the 
edge, and the participation ratio. If the bulk gap were 
the controlling factor then a 0.5 reduction of the band 



gap should shift all observables to smaller disorders, i.e. 
W — >■ W/2. Pane d of Figure [s] shows where the DOS 
would shift to if the bulk gap were the controlling scale. 
We observe that the DOS (and other observables) shift 
slightly in the opposite direction. This confirms that the 
Dirac cone is topologically protected from surface disor- 
der by the bulk band width not the bulk band gap, and 
agrees with the well-established fact that bulk disorder 
destroys the topologically protected surf ace ba nd only 
when it reaches the scale of the bulk band.l^^^ 



IV. DISCUSSION 

The resonance which we have observed in the DOS is 
caused by resonances between individual disordered sites 
and the edge states. As an electron moves along the edge 
of the TI it from time to time becomes almost trapped 
at a particular site and dwells there for a while before 
continuing its journey. This trapping is unable to destroy 
the edge state but it does introduce both a phase shift 
and a localized increase in the surface state's probability 
density. This destroys both the edge state's plane wave 
character and the Dirac cone. The trapping naturally 
also decreases the Fermi velocity of the edge states. In 
pane a of Figure [3] we have plotted the profile of the 
Fermi velocity, which can be easily computed because 
both dimensional analysis and the Dirac relation E = 
VFfi\k\ indicate that it is inversely proportional to the 
density of states. At its peak the resonance divides the 
(average) Fermi velocity by 5.8 and as a result multiplies 
the DOS by the same number. 

In synthesis, our results suggest that surface conduc- 
tion always remains controlled by an effective linear Dirac 
Hamiltonian, albeit renormalized by disorder. The states 
in the bulk gap remain very close to the TI's boundary, 
and we have found that the surface states' energy lev- 
els always remain evenly spaced. These facts argue that 
the in-gap states are described well by a surface Hamil- 
tonian. The robust 2Go conductance indicates that the 
surface states never backscatter, which is consistent with 
spin-momentum locking on a Dirac cone. Moreover their 
DOS is generally independent of energy as expected of a 
2-D Dirac cone. We conclude that disorder renormalizes 
the surface states but does not alter their Dirac physics. 
They are well described by a linear Dirac Hamiltonian 
with renormalized parameters. The main effects of sur- 
face disorder are to reduce the Fermi velocity, increase 
the DOS, and first pull the states into the disordered re- 
gion and later (at very strong disorder) expel them into 
the clean bulk. 
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